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ABSTRACT 


The thermal relaxation process for a spatially unifbrm two 
dimensional plasma in a -uniform dc magnetic field is simulated numer- 
ically. Thermal relaxation times are defined in terms of the time 
necessary for the n\imerically ccmputed Boltzmann H-functlon to de- 
crease through a given part of the distance to its minimum value. 


Dependence of relaxation time on -two parameters is studied: number of 

2 

particles per Debye square a:id ratio of gyrofrequency to plasma 

frequency When becomes » where L is the 


linear dimension of the system, it is found that the relaxation time 
varies to a good approximation as * 
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I. IKTROrOCTION 


In a con^anion paper (Part l)> predictions were derived for 

the rate of thermal relaxation of a two dimensional plasma model in a 

dc magnetic field. The principal prediction which is numerically 

testable is that so long as n/cu > the ratio of gyrofrequency to plasma 

2 . “X/2 

frequency, satisfies the inequality (n/u^) » (^(L/2n^)) y where 

is the Debye length and L is the linear dimension of the system, 
then the relaxation time should be proportional to and to 

(n where is the number density of the charged rods. 

It will be seen below that both of these predictions appear 
to be well fulfilled. In Section 2, the method of simulation is 
briefly described. In Section 3, the results are described and 
carpared with theory. 


2. Method of Simulation 

Tbia dynamical evolution of the plasma was simulated by fol- 
lowing the trajectories of a large number of negatively charged rods 
immersed in a positive neutralizing background. Standard partlcle- 
In-cell techniques with area weighting were enployed (Morse 1970). 
The electric field was determined through Poisson's equation by fast 
Fourier transforms on a 52 x 32 or 64 x 64 grid syst^. Periodic 
boundary conditions were used throughout. The number of rods ranged 
from 4,000 to l6,000. 



k 


Periodically a histogram representing the spatially averaged 
distribution function was calculated. The value of the Boltzmann 
H function was determined by summing the values of 2rrv ^ f(v 
‘ 0in Av^^ obtained from the histogram. The optimum value of 

was found empirically, and for a range of Av^^ (sirqply labeled 
AV in the figures) neither too small nor too large, H(t) was found 
not to depend upon av . (see Figure 1 ). 
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3. Results 

Figure 1 shows a typical plot of H(t) vs t for zero magnetic 
field, for the initial distribution: f(v) = 0 for |y| > 2, f(v) 

= (^tt) for |v| < 2. Times are always expressed in units of the 
reciprocal plasma frequency ^ , and velocities are in units of the 
thermal velocity in the final expected Maxwellian f(y) t_* » 
(2TrVtj^) exp(- v /2V^j^). This makes the natural unit of length 
^th^ ^ ^ Debye length. The particles are loaded randomly in 

position space and in a spiral fashion in velocity. AT is the size 
of the time step. AV is the dimension of a cell size in velocity 
space used in computing the stun which approximates Jdyf0«-f. a = k^/m^ 
is the ratio of Debye length to the cell size m^ used in the particle- 
in-cell conputations . Roughly speaking, larger values of (X mean 
that the short range part of the potential is represented more ac- 
curately. This accuracy is desirable in situations in which short- 
range interactions become dcminant in the relaxation process, as is 

the case when fi/m becomes » 1. 

V 

It is seen from Pig. 1, that while there is a substantial 

-1 

change in the relaxation rate as a decreases from 1.2732 to 0.6366, 
there is no appreciable change upon going to a”^ - 0.3l83. It may 
be safely assumed, then, that the interactions with inpact parameters 
;» 0*5^ are primarily responsible for the decrease of H(t). Figure 1 
should also inspire confidence that in the range in which the com- 
putations are performed, the decrease of H(t) is independent of AT 
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and AV, for al 1 practical purposes. The results shown in Hg. 1 are 
in good agreement with Montgomery and Nielson (l970)» 

Figure 2 shows that for finite magnetic field, increasing a 
continues to increase the speed of thermal relaxation. This reflects 
the fact that in a magnetic field, the thermal relaxation is deter- 
mined hy the close encounters (over-lapping gyroradii). Since 
particle -in-cell ccirputations necessarily "soften" the coulomh 
repulsion at distances ^ a cell size , a decrease in cell size is 
equivalent to a "hardening" of the coulomb repulsion at short dis- 
tances, and thus results in an increased relaxation rate. It is this 
dependence upon the (inaccurately represented) short range part of the 
potential that has, for exaji5)le, precluded our achieving better than 
a factor of two in absolute agreement with numerical calculations of 
energy exchange times for "test" particles and the analytical 
sions for them which can be derived frcm the kinetic equation. 

Figure 5 shows one of otir two main results, the variation of 
the relaxation tine with The relaxation tine ia defined bjr ^ 

setting H = -2.70 (its initial value is -2.551 and i'ts t = ® value is 
-2.857). The = 0 result is close to the value obtained previously 
(Montgomery and Nielson, 1970)* Above abcxit ~ 'we may assert 

that » Wl»/2n^-p)3”^'^^ (our system is l6 units of in linear 

dimension), so that the linear prediction of Part I is appropriate, 
and is seen to be well fulfilled, being accuirately fit by dx^/dn 

= 1.5 for the lower curve and dT./dSl ~ 2.9 for the upper curve. 

KiiiJb 
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Note that the coefficients stand in about the same ratios as 
Between = 0 and =4, we have no satisfactory theory. 

Figure 4 shows the relaxation time r vs the square root of 

RCili 

the nianber of particles per Debye square for two different magnetic 
field strengths. Both sets of data points are seen to be well fit by 
straight lines. 


4. Discussion 

The predictions that the thermal relaxation times of a two 
dimensional magnetized plasma model vary proportionally to D/ and 
have been verified by simulation methods. 
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FIGUEE Ciff>TrOKS 

Figure 1 H-value as a function of time varying (time step), 
(velocity interval), and Ot. The magnetic field 
is zero. 


SYMBOL 

AT 

AV 


X 

0.25 

1/55 

1.2752 

0 

0.50 

1/44 

1.2732 

t 

0.50 

1/55 

1.2752 

+ 

0.25 

1/55 

0.6366 

a 

r 0.50 

1/33 

0.6366 


1 0.50 

1/44 

0.6366 

A 

0.25 

1/55 

0.5185 

A 

1.0 

1/55 

0.3183 


0.50 

1/44 

0.5185 


Figure 2 H-value as a function of time, 
il/ <1^ “ 0. 5 • 

Figure 3 Relaxation time vs magnetic field strength for 
constant 

Figure 4 Relaxation time vs (n^ for constant magnetic 


field strength. 
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